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One dimensional Rectangular, Barlett , Cosine, Hann , 
Hamming, Blackman and Kaiser window are extended to two 
dimensions using circular symmetry and the Hankie transform 
in integral form. It is shown that this transform cannot 
generally be expressed in closed form. Comparisons of the 
above two dimensional windows are then made on the basis 
of two dimensional extensions of the criteria of equivalent 
noise bandwidth, coherent gain, and the so called trade-off 
parameters . 
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I . INTRODUCTION 



Windowing is essential to signal processing and to digi- 
tal filter design. For example, the weighting functions of 
the window are used in harmonic analysis of signals [1] and 
in design techniques for linear-phase digital FIR filter [2]. 

In fact, if an infinite sequence of numbers is truncated, 
then its Fourier transform will be subject to the well-known 
Gibb's phenomenon; which will cause over shoot and ripples 
before and after any point of discontinuity in the frequency 
domain. The overshoot and ripples will not disappear no 
matter how long the sequence is . 

To overcome this undesirable effect in the frequency do- 
main, the finite sequence should be modified using a weighting 
function, called a window. For example, if w(n) is a weight- 
ing sequence and h(n) is an infinite sequence of number, then 
the windowed finite sequence is: 

h (n) = w(n) . h(n) (1.1) 

w 

where : 

f(n) |n|<N 2) 

^ ^ 0 elsewhere 

In this example the weighting functions performs both 
truncation and modification of the signal at the same time. 
Multiplication of two sequences in time domain is equivalent 
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to the convolution of their transforms in the frequency 
domain. Consequently the windowing method has the effect 
of smoothing out the rough points - moderating the overshoot 
and ripples - in the original frequency response. 

The technique of windowing is shown graphically in 
Fig . 1.1. 

There are many window functions available in many pieces 
of literature. [1, 3] According to the way in which the 
window function is formed, they can be separated into two 
categories ; basic windows and constructed windows . A con- 
structed window is formed from a product, a sum, or a con- 
volution of two or more basic windows. 

Seven best-known basic windows are listed in Table 1.1, 
and they are: Rectangular, Barlett, Cosine, Hann, Hamming, 

Blackman and Kaiser window. Most of the windows are named 
after individuals. These seven basic windows will be dis- 
cussed thoroughly in later chapters. 

Although Huang [4] shows that a good two dimensional 
window can be derived from one dimensional window, most 
windows found in the literature are given in only one-di- 
mensional form. It is the purpose of this thesis to present 
and discuss two-dimensional windows derived from those one 
dimensional form, using Huang’s theorem. 

The easiest way to compare various windows is to create 
some appropriate figures of merit for windows. This thesis 
presents those figures of merit that are used to compare the 
given window functions. 
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A. TECHNIQUE OF WINDOWING [2] 



A graphical interpretation of windowing is given in 
Fig. 1.1 (a non-recursive filter design using the window 
method) . 

The first figure shows an infinite impulse response of 
a digital filter, h(n) , and its periodic frequency response 
H(w) . 

The next figure shows a rectangular window r(n), with 
its frequency response, R(w) . 

The third figure consists of multiplication of h(n) and 
r(n), which produces a truncation of the filter impulse 
response, h^(n). Multiplication in the time domain is equi- 
valent to convolution in the frequency domain. The convolu- 
tion of H(w) and R(w) is H^(w) containing overshoot and 
ripples in its frequency response of the corresponding digi- 
tal filter as shown in the figure. 

Truncation of impulse response is required for a non- 
recursive digital filter design, since the filter produces 
an output which depends only upon the present and a finite 
number of past inputs . 

The fourth figure shows the non-rectangular window w(n) , 
and its periodic Fourier Transform W(w). The purpose of 
this window is to alter the truncated filters impulse response, 
h^(n), slightly in order to minimize the overshoot and ripples 
in its frequency response H^(w) . 

The last figure shows the finite impulse response h (n) , 

w 

which consists of multiplication of h^(n) and w(n). The 
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frequency response of h (n) is produced by convolving H (w) 

with W(w), resulting in a smooth frequency response H (w) . 

w 

The technique of windowing applied to a filter design, 
results in an improved filter with a smooth frequency response. 

B. THESIS ORGANIZATION 

The following discussion presents the organization of 
this thesis. 

As mentioned before, the best way to discuss windows is 
to define parameters which can be used to characterize a 
window. Logically, the window figures of merit are first 
discussed in the next chapter. Chapter II. 

Chapter III gives an overview of one-dimensional windows 
and specifies their characteristics, using parameters defined 
in Chapter II. Table III.l, given at the end of this chapter, 
tabulates one-dimensional windows and their parameters. 

The transformation from one-dimensional to two-dimensional 
is given in Chapter IV, where we present three major discus- 
sions. First we present the expansion of the Fourier Trans- 
form from one-dimension to two-dimensions. Then we present 
how Huang [4] expanded a good one-dimensional window into a 
good two-dimensional window. Finally, the definition of one- 
dimensional figures of merit are expanded into two-dimensional 
window figures of merit , 

Chapter V presents two-dimensional windows in detail, and 
specifies their parameters as defined in Chapter IV. Table 
V.l, given at the end of this chapter, tabulates two-dimen- 
sional windows and their parameters. 
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h(n) 




r (n) 








1 R(w) 







w 








Figure 1.1 

Application of Window Function To 
FIR Filter Design 
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Two-dimensional Fourier Transform of window functions 
are hard to present in a closed form. Fortunately a closed 
form Fourier Transform is not required in studying windows, 
so two-dimensional Fourier Transform of the windows are 
presented in integral form. The effort to present a two- 
dimensional Fourier Transform of a window in a closed form 
is given in the Appendix. 
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TABLE I.l 



Seven Basic One-Dimensional Window Functions 



No. 


Window 


Function 








1 


Rectangular 


W(x) = 1.0 


+ M 

a 


|xl 


1^ 


2 


Barlett 


O 
1 — 1 

II 

X 


|x| 


1 <a 


3 


Cosine 


W(x) = cos 


X 

2a 


|x| 


1 <a 


4 


Hann 


W(x) = I + 


1 

-T- COS — X 

2 a 


|x| 




5 


Hamming 


W(x) = 0.54 


+ 0.46 cos — X 
a 


|x| 


|j<a 



Blackman 



W(x) = 0.42+0.5 cos ^ X + 0.08 cos — x 

2 a 



X <a 



Kaiser 



W(x) 



I (S / 






X <a 
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II. WINDOW FIGURES OF MERIT [1, 3] 



Since windows are used in frequency response analysis to 
reduce undesirable effects - overshoots and ripples in the 
frequency domain - it is very important to find the para- 
meters or figures of merit of the windows that can be used 
to identify their characteristics and to indicate their per- 
formance . 

There are many convenient measures that can be used as 
parameters . 

The common desirable characteristics of a window are: 

1. A narrow main lobe width of the frequency response 
of the window, containing as much of the total energy as 
possible . 

2. Side lobes of the frequency response decreasing in 
energy rapidly as frequency increases. 

These two common characteristics certainly can be used as 
parameters to identify the window's performance. 

Harris [1] and Geckinli [3] define one-dimensional window 
performance parameters which will be repeated here. These 
parameters are used throughout this thesis and will be expan- 
ded for the two-dimensional case in later chapters. The 
parameters are; equivalent noise BW, coherent gain, and 
trade-off parameters. 
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A. EQUIVALENT NOISE BANDWIDTH (ENBW) 

The ENBW parameter is the width of a rectangular filter 
in the frequency domain with the same amplitude as the peak 
power gain of the window, and with a width such that the 
total noise power in the rectangular filter is equal to the 
total window noise power. 

A graphical explanation of the parameter is shown in 
Fig. 2.1. The first figure shows a window frequency response 
which consists of a main lobe and side lobes. The second 
figure shows a discrete frequency component of a signal, 
together with white noise. The third figure shows the result 
of windowing . 

The window's frequency response convolves with the dis- 
crete component of the signal plus noise. In this case the 
window acts as a filter, gathering information over its 
bandwidth. The last figure represents the area of the windowed 
discrete frequency response of the signal plus noise, at dis- 
crete frequency W^. 

It is clear from the figure that to minimize the accumu- 
lated noise collected by the window, it is required that the 
window has a narrow bandwidth. 

A convenient measure of window's bandwidth is called 
equivalent noise bandwidth. 

Noise power accumulated by the window is defined as: 

N 

Noise Power = ^ /“ |W(w) r dw (2.1) 

where N^ = noise power per unit BW. 
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W W 

o 



Figure 2.1 

Graphical Interpretation Of Convolution Of A Window 
And Spectral Amplitude Of A Signal, At w=w^ 
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M II 




Parseval's theorem allows equation (2.1) to be computed as: 



Noise Power = 2^ W^(t) dt (2.2) 

The integral in (2.1) and (2.2) represents the area under 
the window curve . If the area is divided by height , then 
the result is the width of a rectangle that has the same area 
with the window. 

The peak gain of signal occurs at w=0, and is defined as: 

T 

Peak signal gain = w(0) = / w(t) dt (2.3) 

-T 

T 

Peak power gain =w^(0) = [ /w(t) dt]^ (2.4) 

-T 



Then the ENBW is equal to the noise-power divided by peak 
power: 



T 

/ w^(t) dt 

ENBW = ^ (2.5) 

[ / w(t) dt] 

-T 



The above expression is a normalized ENBW, after it is 

divided by factor ^ . This parameter will be used as a 

2T 

window figure of merit. 



B . COHERENT GAIN 

Coherent gain is best defined using a sinusoidal signal. 
In Fig. 2.1 a component of signal frequency response at a 
discrete frequency is shown. Convolving this signal with a 
window’s frequency response can be thought as passing the 
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W(w) 



2 




Figure 2.2 

Equivalent Noise BW Of A Window 
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signal through a filter. Using this view, the coherent gain 
of a window can be examined by passing a sinusoidal signal 
through a filter. 

Let the sinusoidal signal be defined as: 

f(x) = Ae^”^ + n(x) (2.6) 

where n(x) = white noise. 

The filter output on the windowed signal is given as: 

f(x)=Aw(x)e^^^ + n(x)w(x) - a<x<a (2.7) 

w 

The frequency response of the windowed signal only (less noise) 
is given as: 

F (f) = /A W(x)e^"^ e-^"^ dx = 

^ -a . 

a 

= A / W(x) dx (2.8) 

-a 

The windowed frequency response is proportional to the input 
amplitude A. The proportionality factor is the integral of 
the window function. 

For a rectangular window this factor is 2a. For any other 
window the factor is reduced due to the window smoothly going 
to zero near the boundaries. 

The coherent gain is defined as the proportinality factor 
divided by factor 2a, the extension of the window (normalized). 

COH.GAIN = ^ / W(x) dx (2.9) 

L d. 

-a 
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The coherent gain is the second parameter that will be used 
as window figure of merit. 



C. TRADE-OFF PARAMETERS 

Trade-off parameters are best explained using Fig. 2.3. 
The top figure shows the Hamming window in the time or spa- 
tial domain, and the bottom figure shows its Fourier trans- 
form, which is given in log frequency versus normalized log 
magnitude (20 log W(f)/W(0)). 

The window's characteristics in frequency domain are 
largely determined using these trade-off parameters. 

From Fig. 2.3(b) these parameters are: 



1 . a^ : 


highest side lobe level in dB. 


2 . a^ : 


side lobe level at a higher frequency, 
where asymptotic assumption is valid for 
the window behavior. 


3 . b : 


frequency at which main lobe drops to the 
peak side lobe. 


4. d : 


side lobe fall-off rate in dB/octave. 



Clearly from the above definitions, the smaller the value 
for b, the better the window performance*, the smaller the 
value for a^, the smaller the leakage through the near side 
lobes ; the smaller the value for a .^ , and the larger the value 
of d, the smaller the leadage through the far side lobes. 
These parameters are expanded for two-dimensional case in 
a later chapter. 
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Plot Of Hamming Window, Is Used To Illustrate The 
Definition Of Trade-Off Parameters a^ , b and d. 
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III. ONE DIMENTIONAL WINDOW FUNCTIONS 
[1, 3, 5, 6, 7] 

The seven common window functions listed in Table I.l, 
are reviewed in detail in this chapter. Figures of merit 
for each window are presented using the definitions given in 
the preceding chapter. The window functions treated here are 
in continuous form rather than in discrete form. 

Without loss of generality, let the window function, 
w(x) , be unity at the origin and limited to the interval 
|x|_< a. Presented in mathematical form this means: 

00 

w(o) = . f W(f) df = 1 (3.1) 

— 00 

w(x) = 0 for |x|^a (3.2) 

The Fourier Transform of the above window is given as: 

a 

W(f) = / w(x)e dx |x|_<a (3.3) 

-a 

The transforms of the windows generally consists of a main 
lobe with a large peak amplitude and side lobe with compara- 
tively small peak amplitudes. To better display the presence 
of the side lobes, logarithmis scales are used for plotting 
in place of linear scales . 
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A. RECTANGULAR WINDOW [1] 

The rectangular window is commonly called a boxcar or 
Dirichlet window. Its main purpose is to truncate the 
length of a function to a finite extent, without affecting 
its amplitude. 

The expression for this window is given as; 

w(x)=1.0 (3.4) 

The Fourier Transform of the window is given as: 

W(f) = i dx = 

Z 0. 

-a 



sinaw 

aw 



sine (aw) 



(3.5) 



The transform is seen to be a Dirichlet kernel. In this 
report the extent of the window is taken as a=0.5, and the 
ENBW is given as : 



0.5 

f dx 

ENBW = ^ = 1.0 . 

[ / dx ] 

-0.5 



The coherent gain is given as: 

0.5 

COH.GAIN = / dx = 1.0 

-0.5 
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Plots of the window function and its tranform are given 
in Fig. 3.1(a) and (b) respectively. The main lobe width 
(the width between origin and the first zero crossing) is 
1.0, and the first side lobe peak is approximately 13 dB 
down from the main lobe peak. The side lobe fall off is 
6.0 dB/octave. 

B. BARTLETT WINDOW 

The Bartlett window is also called Fejer or triangular 
window, and can be thought of as a convolution of rectangu- 
lar window with itself. This window is the simpliest of all 
windows that changes the data values by multiplication by a 
number other than unity. The mathematical expression of this 
window is given as: 



w(x) = 1.0 - 

3 . 




(3.6) 



The transform of this window is given as: 



W(f) = - 

3 . 



n I ^ U “ jwx , 

/ (1.0 - — ^)e dx 

3 . 



a 



1 



o 



/ (1.0 + dx -H 

a 



a 



1 




a 



1 




a 



/ X e dx + 



a -a 




2a o 
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- 0 . 5 



000 



0 . 5 



X 



(a) 




(b) 

Figure 3.1 

(a) Rectangular Window, (b) Spectral Window 
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T a , - ] wx 

1 . de 

- — o / X : 

2 -1 w 

a o 



n law . -iwx o 



a -]w 



~~2 ' -iw 
a -a 



- / 



° e-]™dX) ^ 



-a 



-3W a -]w 



T -1WX a 

7 - 3 “ o 



a - jwx 
/ - — = dx } 



, law - 1 aw , law - i aw 
l,e-^ -e V 1 /e-’ -e 






., . ± c 



T O -1WX 

1 , e 



]W 



) - - 5 - / 

a -a 



- jw 



dx + 



a -]wx 

+ — 75- / 1 dx 

2 -1 w 

a o 



1 T 2 . o , T 2 . 

1 /lx -]WXr, 1 /lx ^-]WX 

- ( : ) e + —ry ( r— ) e 

2 -iw ' 2 -iw 

a a a 



+ 1 (1 . + 1) = -4^ (e3^“ - 2 -f 

aw aw 



T 2 jaw/ 2 -jaw/ 2 

"fj > 



.sin aw/ 2 x , . / o \ 2 



( 3 . 7 ) 



The transform is seen to be the squared of Dirichlet Kernel. 
If the extent of the window is taken as a= 0 . 5 , then the 
parameters ENBW and COH.GAIN are: 
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ENBW 



0.5 2 

f (1- 7 t\) dx 

-0.5 



0.33 
2 “0,25 



1.33 



0,5 

[ / (1 

-0.5 




0.5 

COH.GAIN = / (1 

-0.5 



dx = 0.5 



The window function and its transform are shown in Fig. 
3.2(a) and (b) respectively. Its main lobe width is 2.0, 
twice that of the rectangular window, and the side lobe fall- 
off is 12 dB/octave. 

C. COSIME WINDOW 

Cosine window is actually one member of a family window 
with a general form given as: 



Type of this window depends upon parameter a, which usually 
has a value between 1 and 4. Of particular interest in this 
family window is Hann window, where a=2, which will be con- 
sidered later. 

The cosine window is defined with a=l, as: 



wCx) = cos x) 




(3.8) 



w(x) 





(3.9) 



The Fourier Transform of this window is given as: 
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H ^w(x) 



O 





Figure 3.2 

(b) Spectral Window 



(a) Barlett Window, 



1 3. . , TT X 

2a ^ 

-a 



2 a dx + -5— / e 

Z d. 

-a 



1 ^ -~j(w+-;^)x 



2a dx 



, sin a(w— sin a(w+-^) 

1 ,2a . 1 2a 

2 / TT s 2 _ / ^ . IT 



a(w— 



a(2.^) 



1 • / TTv l • /j.TT\ 

TT sine a(w — ^ sine aCw+^r-) 

2 2a 2 2a 



(3.10) 



The transform of this window is a summation of two shifted 
Diriehlet Kernel. 

If the extend of the window is taken as a=0.5, then the 
parameters ENBW and COH.GAIN are given as: 



ENBW 



u . o 2 

/ (eosTTx) dx 
-0.5 

0.5 2 

C / eosTTx dx] 

-0.5 



0.5 

0.404 



1.23 



COH.GAIN = / eosTTx dx = — = 0.64 

-0.5 ^ 

Plots of the window and its transform are given in Fig. 
3.3(a) and (b) respeetively . From the transform plot, the 
main lobe width is 1.5 times that of the reetangular window. 
The first side lobe peak is approximately 23 dB down from 
the main lobe peak, while the side lobe fall off rate is 
12 dB/octave. 
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w(x) 



O : 




0 log „ |W(f )/W(0) I 

b=l. 35 20. log f 




Figure 3.3 

(a) Cosine Window, (b) Spectral Window 



D. HANN WINDOW 

The Hann window is one member of the family windows given 
ct 

by cos_x, where a=2. The window equation is written as: 

w(x) = cos = 2*2 "aT" (3,11) 



The Fourier Transform of this window is given as 



W(f ) 



i / 

2a ^ 

-a 






TT 

cos— x) 



-1WX 



dx 



, a . 1 ^ 

= T . — / e dx + r — f cos —X e dx = 

4a 4a a 

“■a ““a 



1 sin aw 

2 - aw 



+ 



1 

4 



, TT . 

Sin a(w — ) 
^ 

a( w — ) 
a 



+ 



1 

4 



sin a(2+— ) 
a 

a(w+— ) 
a 



(3.12) 



The Transform of this window consists of a summation of 
three Dirichlet Kernels, one at the origin and the two others 
shifted — radians from the origin. 

The three terms on the transform of this window are 
shown in Fig. 3.4, where the shifted Dirichlet Kernels are 
located at the first zeroes of the center kernel. The peak 
amplitudes of the shifted kernels are half that of the center 
kernel; also the side lobes of the translated kernels are 
about half the size of, and of opposite phase from those of 
the center kernel. The summation of the three side lobes in 
phase opposition, tends to cancel the overall side lobe structure. 
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Figure 3.4 

Transform Of Hann Window As A Sum 
Of Three Dirichlet Kernel 
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The cancellation of the side lobes structure suggests 
a constructive technique to define new windows . The best- 
known windows with this type of side lobes cancellation are 
Hamming and Blackman windows , which are given in the next 
section . 

If the extent of the Hann window is taken as a=0.5, the 
parameters ENBW and COH.GAIN are given by: 



/ (h+h coSfJ^)^ dx 
ENBW = 

[ / (h + h COS^X) dx]^ 

-0.5 



0.375 

0.25 



1.5 



0 . 5 

COH.GAIN = / (h+h cos ) dx = 0.5 

-0.5 

The window and its Fourier Transform are plotted in Fig. 
3.5. The main lobe width is 2.0, twice that of the rectan- 
gular window and the first side lobe level is 23 dB down 
from the main lobe. The side lobe level fall off rate is 
18 dB/octave. 



35 




I 




rH W(X) 
o . 





Figure 3 . 5 

(a) Hann Window, (b) Spectral Window 
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E. HAMMING WINDOW [5] 

Hamming window can be thought an an improvement of Hann 
window, in a sense that the first side lobe cancellation 
is even greater. From Fig. 3.4, the summation of the three 
side lobes produce an inexact cancellation. The cancellation 
can be made more pronounced by adjusting relative size of 
each term. 

Exact cancellation of the side lobes was first shown by 
Hamming using a window given in general form as: 

w(x) = a + (1-a) cos — x |x|<a (3.13) 

where a is a coefficient that will be defined later. The 
Fourier transform of the. window is given as: 

W(f) = asinc aw + 0.5(l-a)sinc a(w-— ) 

3 

+ 0.5(l-a)sinc a(w+— ) (3.14) 

3 

From (3.14) the best cancellation of the first side lobe 

25 25 

occurs when a = yg- , (-jg- = 0.543478261). If a is selected as 

25 

0.54, an approximation of then an improvement in cancel- 

lation of the side lobes occurs. The value a=0.54, gives 
rise to a window called the Hamming window, which is written 
as : 



w(x) 



0.54 + 0.46 cos 



TTX 



|x|^ 



(3.15) 
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The Fourier transform then is given as: 



W(f) = 0.54 sine aw + 0.23 sine a(w-— ) + 

3 . 



+ 0.23 sine a(w+— ) 

3 



(3.16) 



If the extent of the window is a = 0.5, then the para- 
meters ENBW and COH.GAIN are given as: 



ENBW = 



0.5 IT 2 

/ ( 0 . 54 + 0 . 46eos — rx) dx 

-0.5 



0.3974 _ 



0.5 



[ / ( 0 . 54 + 0 . 46eoSpT^x) dx] 

-0.5 



0.2916 



= 1.36 



0.5 2 

/ ( 0 . 54 + 0 . 46oos^y^x) dx 

-06 ^ ^ 03974 

COH.GAIN = ^ = 1.36 

[ / ( 0 . 54 + 0 . 46eos-pr^x) dx] 

-0.5 

The Hamming window and its transform are ploted in Fig. 
3.6(a) and (b) respeetively . From its transform plot, the 
main lobe width is twiee that of the reetangular window, and 
the first side lobe level is 43 dB down from the main lobe. 
The side lobe fall of rate is 6 dB/octave. 

F. BLACKMAN WINDOW [6] 

Hann and Hamming windows are examples of windows con- 
structed as the summation of shifted sinc(x) functions. 

Blackman suggested a general rule by which window 
functions could be constructed as follows : 
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Figure 3 . 6 

(a) Hamming Window, (b) Spectral Window 



39 



r 




k 

w ( X ) = Z 

n ni 



cos m ( — ) 
a 



I ^ I 



(3.17) 



and subject to the constraint; 



k 

Z a = 1.0 
_ o m 



(3.18) 



The Fourier transform of this window becomes: 



k 



a 



W(f) = Z (-1)^ {sine a(w-— m) + sine a(w+— m)} (3.19) 

m=0 ^ ^ 



being non zero. Their spectral windows are summations of 
three shifted kernels. 

Any value of k will produce a new window with a summation 
of (2k-l) kernels in its tranform, and k non zero corficients, 
however one way to achieve windows with a narrov; main lobe 
is to restrict k to a small interger. 

Blackman examined this window with k=3, and found best 
cancellation of the third and fourth side lobes occurs when: 



Hann and Hamming windows are of this form with a^ , and a^^ 




9240 

^1 ■ 18608 



0.496 560 62 = 0.50 



1430 _ 

^2 ■ 18608 



0.076 848 67 = 0.08 
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The window function which uses these three coeficients 
is known as the Blackman window, and is given as: 



w(x) = 0.42 + 0.5 cos— + 0.08 cos— x I x I <a (3.20) 

a a ' ' — 



The Fourier transform of this window is: 



W(f) = 0.42 sine aw + 0.25 sine a(w-— ) + 0.25 sine a(w+— ) 

a a 

9 tT 9 7T 

+ 0.04 sine a(w ) + 0.04 sine a(w+— ) 

a a 

(3.21) 



If the extent of this window is a=0.5, then the ENBW and 
COH.GAIN parameters are given by: 



ENBW = 



0.5 7 2 

/ (0.42 + 0.5 cos —X + 0.08 cos — ^x) dx 

-0.5 ® a 

“•5 n 2ir 2 

[ f (0.42 +0.5 cos— X + 0.08 cos — x) dx] 
-0.5 ^ ^ 



0 . 2646 
0.1764 



= 1.57 



COH.GAIN = 



0.5 . 

f (0.42 + 0.5 cos — s + 0.08 cos — x) dx = 0.46 
-0.5 ^ ® 



The window and its transform are shown in Fig. 3.7(a) and 
(b) respectively. The main lobe width is 3.0, three times that 
of the rectangular window, and its first side lobe level is 
41 dB down from the main lobe. The side lobe fall off rate is 
18 dB/octave. 
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G. KAISER WINDOW [7] 

The Kaiser window is also called a Kaiser-Bessel window 
and its mathematical expression is given as: 



w(k) 



I [3 /i.o-(-)^] 

O 3. 



|x|_^ (3.22) 



where I (x) = zeroth order modified Bessel-f unction , and is 
o ’ 

defined as: 



I (x) 
o 



1 + 



S C 

n=l 



^2’ 



1 

n 



2 

] 



(3.23) 



The Fourier transform of this window is given as: 



W(f) = 



2 

I (3 ) 



sinh[3 






f <fa 



(3.24) 






where fa = 3 / tt 



From (3.24) if f=fa, the term inside the square root 

2 

becomes zero and W(fa) = j f>fa, the term becomes 



i /(i_)2 1 
^ / ^fa^ 

becomes : 



and using the identity sinh ix = i sin x, (3.24) 



sin 



W(f) = 



C3 / (^)^-l] 



■fa^ 



I„(B) 



f >fa 



(3.25) 






' fa 



) -1 
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The usual range of 3 is between 4 and 9 , corresponding to 
a range of relative peak heights from 3.1% down to 0.047%. 

As 3 increases, the side lobes decrease in amplitude while 
the main lobe widths increase. By varying 3 the proportion- 
ality between the width of the main lobe and the height of 
the side lobe can be adjusted as desired. 

When it is given 3=6 and a=0.5 then the parameters ENBW 
and COH.GAIN are given as: 



ENBW 



0.5 

/ 

- 0 . 5 



0 . 5 

C / 

- 0 . 5 



I^C6 

/ O U • b V 

^ ^ 

O 

O U . 0 

TTeT 

o 




dx]^ 



0.36674 
" 0.25028 



1.5 



COH.GAIN 



0 . 5 
/ 

- 0 . 5 



dx 

I [6] 
o 



0.5 



Fig. 3.8(a) and (b) shows the window and its transform 
respectively. The main lobe width is 2.0, twice that of 
the rectangular window, and the first side lobe level is 
43 dB down from the main lobe. The side lobe fall off rate 
is about 6 dB/octave, 

For 3=6 and 3=9 the window resembles to Hamming and 
Blackman window respectively. 
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Figure 3.8 

(a) Kaiser Window, (b) Spectral Window 
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TABLE III.l 



ONE- 


-DIMENSIONAL 


WINDOW AND 


FIGURES OF 


MERIT 




WINDOW 


ENBW 


COH.GAIN 


^1 


^2 


b 


d 








dB 


dB 


Hz 


dB/ octave 


Rectangular 


1.0 


1.0 


-13 


-35 


0.81 


-6 


Barlett 


1.33 


0.5 


-26 


-58 


1.63 


-12 


Cosine 


1.23 


0 . 64 


-23 


-65 


1.35 


-12 


Hann 


1.5 


0.5 


-32 


-88 


1.87 


-18 


Hamming 


1.36 


0 . 54 


-43 


-53 


1.91 


-6 


Blackman 


1.57 


0.46 


-58 


-94 


2.82 


-18 


Kaiser 


1.5 


0.49 


-43 


-66 


2 . 2 


-6 
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IV. TRANSFORMATION FROM ONE-DIMENSION 
TO TWO-DIMENSION [4, 8, 9] 

The two-dimensional Fourier transform, two-dimensional 
window functions, and window figures of merit are discussed 
in this chapter. 

A. TWO-DIMENSIONAL FOURIER TRANSFORM [8] 

The Fourier transform of f(x,y), a function of two inde- 
pendent variables x and y, is .represented by F.{f(x,y)} or 
F(u,v), as is defined as: 

00 00 

F{f(x,y)} = F(u,v) = / / f (x,y)e"^^'^^'^^y^dxdy (4.1) 



The transform is itself a function of two independent 
variables u and v, which are generally referred to as fre- 
quencies . 

Similarly the inverse Fourier transform of a function 
F(u,v) is represented by F ^{F(u,v)} or f(x,y) and is 
defined as: 



00 00 

F~^{F(u,v)} = f(x,y) = / f F(u,v)e^ dudv (4 



2 ) 



As mathematical operations, the transform and its inverse are 
very similar differing only in the sign of the exponent 
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appearing in the integral. 

In order that the integral exist, the constraints on 
the function f(x,y) are given as: 

1. f(x,y) must be absolutely integrable over the infinite 
x,y plane. 

2. f(x,y) must have only a finite numbers of discon- 
tinuities and a finite number of maxima and minima 
in any finite extent. 

3. f(x,y) must have no inf inite discontinuites . 

The basic definition of the Fourier transform pair leads 
to a rich mathematical structure associated with the trans- 
form operation. Considered below, a few of the basic mathe- 
matical properties that will find use in later discussion. 
These properties are presented as mathematical theorems, 
followed by a brief statement of their physical significance. 

1. Linearity theorem: 

F {af+$g} = a^{f}+B^{g} (4.3) 

The transform of a sum of two functions is simply the sum 
of their individual transforms, 

2. Sealing theorem: 

F {g(ax,by)}= G (^, (4.4) 

A stretching of the coordinates in the space domain (x,y) 
results in a construction of coordinates in the frequency 
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domain (u,v), plus a change in the overall amplitude of the 
spectrum. 

3. Shifting theorem: 

F {f(x-a, y-b)} = F(u,v)e”^ (4.5) 

A translation of a function in the space domain intro- 
duces a linear phase-shift in the frequency domain. 

4. Parseval's theorem: 

00 00 OO 00 

/ / |f(x,y) dxdy = / / |F(u,v)| dudv (4.6) 

_00 _ 0 O ->00 —00 



This theorem is interpretable as a statement of conserva- 
tion of energy . 

5. Convolution theorem: 

00 OO 

F { / f f ( z ,n)g(x-z ,y-n)dzdn} = F(u,v) G(u,v) (4.7) 

— OO —00 

The convolution of two functions in the space domain is 
entirely equivalent to the simpler operation of multiplying 
their individual transforms. 

6. Autocorrelation theorem: 

00 00 o 

F { / / f(z,n)f“(x-z,y-n)dzdn} [F(u,v)r (^-8) 
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similarly : 



00 00 

F {|f(x,y)| } = / f F(z,n)F ( z+u ,n+v) dudv (4.9) 

..00 ~00 

These theorems are special cases of the convolution 
theorem. 

7. Fourier identity: 

At each point of continuity of f(x,y): 

F ^ xe {f(x,y)} = F F ^{f(x,y)} = f(x,y) (4.10) 

The successfive transformation and inverse transforma- 
tion of a function yields that function again, except at 
point of discontinuity. 

1. Separable and Circularly Symmetric Functions [9] 

A function of two independent variable is called 
separable with respect to a specific coordinate system if it 
can be written as a product of two functions. 

A function f(x,y) is separable in the rectangular 
coordinate x,y if: 

f(x,y) = f^(x) f2 (y) (4,11) 

A function f(r,9’) is separable in the polar coordinate 

r , & if : 

f(r,&) = f^(r) f^(^) (4.12) 
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Separable functions are more convenient to deal with 
than more general functions . Separability allows complicated 
two dimensional manipulations to be reduced to simpler one 
dimensional manipulations. 

A function separable in rectangular coordinates has 
the simple property that its Fourier transform is obtained as 
a product of two one dimensional transforms 



00 

F {f(x,y)} = / 

— 00 



" f(x,y)e-3'=‘“*''>'>dxdy 



00 

= f 
— 00 



^00 

f ^(x)e”^'^^dx / f2(y)e”^^^dy 

— 00 



= F (f^(x)} F {f 2 (y)l (4.13) 

Thus the transform of f(x,y) is itself separable into 
a product of two factors. The process simplifies to a familiar 
one dimension manipulation with a function of u, and a function 
of V . 

Functions separable in polar coordinates are not so 
easily handled as those separable in rectangular coordinates. 
But it is generally possible to demonstrate that two dimen- 
sional manipulations can be performed by means of a series of 
one dimensional manipulations. The simplest class of functions 
separable in polar coordinates are those functions which ex- 
hibit circular symetry. 
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A function is said to be a circularly symmetric if it 
can be written as a function of radius r, independent of the 
variable 9 . 

f(r,9) = f(r) (4.14) 

Such functions are very important, since a window function can 
be extended into a two dimensions if it is given a circularly 
symmetric form. Accordingly, special attention should be 
paid to these functions. 

The Fourier transform of a function in rectangular 
coordinates is given by: 



F(u 


00 00 
,v) = f f 

~oo »oo 


f(x,y)e ^ ^'^^"'"'^^^dxdy 


(4.15) 




If both x,y 


and u,v coordinates are 


transformed into 


a polar 


coordinate 


system then: 




r = 


/ 2^ 2 
/ X +y 


x=rcos9 




9 = 


tan (— ) 

y 


y=rsin9 




and 

P = 


y 2^ 2 

✓ u +v 


u=pcos4> 




4> = 


tan (-) 


v=psin<j) 





Substituting the polar coordinates into (4.15) and 
using the transform identity: 
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it; 



h 



dxdy = rdSdr 



ux+vy = pcos(j)r cos4> + psin(j)r sin0 = r cos(0~4)) 



The Fourier transform in (4.15) becomes: 



F(p,4>) = / 
0 



00 



2tt 

f r f (r , 0 )e 
0 



-jpr cos(e-4.) 



(4.16) 



If the Bessel function identity: 



J ( pr ) 
o 




(4.17) 



is used to simplify (4.16), then the transform becomes: 



function is itself circularly symmetric, and can be found by 
performing a one dimensional manipulation as per (4.18). The 
expression in (4.18) is known as the Fourier-Bessel transform 
or as the Hankel transform. 

By identical arguments, the inverse Fourier transform 
of a circularly symmetric function F(p) can be expressed as: 



00 



F(p) = 2 tt / r f(r) J (pr) dr 

n o 



(4.18) 



Thus the Fourier transform of a circularly symmetric 



00 



f(r) = 2 tt / pF(p) J (rp) dp 
0 ° 



(4.19) 
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There is no difference between the transform and its inverse 
operation in the case of Hankel transform. Using the notation 
H { } to represent the Hankel transform operation: 

H {f(r)} = H H {f(r)} = f(r) (4.20) 

From the Hankel transform pair (4.18) and (4.19), two 
dimensional Fourier transform may be found from one dimen- 
sional mainpulations . 

Huang [4] has shown that a good two dimensional window 
function can be derived from a good one dimensional function, 
as long as it is circularly symmetric. This is discussed in 
the next section. 

The Hankel transform in (4.18) will be used to cal- 
culate the Fourier transform of a two dimensional window. 

B. TWO DIMENSIONAL WINDOW FUNCTION [4] 

If W 2 (u,v) and vi^(.x,y) is a two dimensional Fourier 
transform pair, then its inverse transform relation is given 
as : 



w^(x,y) 
For y=0 



W2(x,o) 



7 7 W2(u,v)e3‘'^^*''y>dvdu 

—oo — oo 

(4.21) becomes: 

OO oo 

f { f W 2 (u,v) dv}e^^^ du 

— OO — oo 



(4.21) 



(4.22) 
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If W 2 (x,y) is given as a circularly symmetric function, 



7 2 2 

then W 2 (x,y) = w( / x +y ) . For circularly symmetric 


function 


W 2 (x,o) = w(x) 

And its inverse Fourier transform is given as: 


(4.23) 


00 

w(x) = f W(u)e^^^ du 
.00 

From (4.22) and (4.23) it is evident that: 


(4 . 24) 


00 

W(u) = / W 2 (u,v)dv 

— 00 


(4.25) 



The expression in (4.25) is called a profile function, 



and similarly the profile function in v is expressed 


as : 


00 

W(v) = / W 2 (u,v)du 

.00 

Now consider the following two equations: 


(4.26) 


1 for u > o 

H(u) = { 

0 u < o 

and 


(4.27) 


1 for u > o for all v 

H2(u,v) = { 

0 u < o for all v 


(4.28) 
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The expressions in (4.27) and (4.28) are one-dimensional and 
two dimensional step function respectively. 

If W 2 (u,v) is a two dimensional window’s Fourier trans- 
form, then its convolution with the step function is given 
as : 



OO CO 



W^(u,v)9 H_(u,v) = / f W„(t,s)H(u-t , v-s)dtds 



OO OO 



U OO 



/ { / W (t,s)ds}dt 

O -OO 



(4.29) 



where t,s are dummy variables. 

Using the profile-function identity given in (4.25) and 
(4.26), it can be shown that: 



W(t) = / W (t,s)ds 

— 00 



(4.30) 



Equation (4.29) becomes: 



W2(u,v) 9 H2 (u,v) 



u 

/ W(t)dt 
o 



(4.31) 



Expanding the right hand side of (4.31) into a convolution form 



u 



/ W(t)dt = / W(t) H(u-t)dt 



= W(u) 0 H(u) 



(4.32) 



56 



I 





Then substituting (4.32) to (4.31): 



0 H 2 (u,v) = W(u) 0 H(u) 



(4.33) 



Huang states that if a two dimensional function is given 
in circularly symmetric form, then its Fourier transform 
satisfies the relation given in (4.33). 

To show the benefit of the results in (4.33) consider an 
example of designing a low pass filter. 

EXAMPLE : 

The ideal Fourier transform of one dimensional Low Pass 
Filter is given as: 



The impulse response f(x) of the filter is equal to its 
inverse Fourier transform, which has an infinite duration. 

In reality, the impulse response has to be finite. By using 
a window function, the infinite duration of the impulse res- 
ponse is truncated to become a finite impulse response as 
given as : 



1 for |u| _< A 



F(u) = { 



(4.34) 



0 



u > A 



g(x) = f(x) • w(x) 



(4.35) 



where : 



w(x) = 0 for |x|>'a 



(4.36) 
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The windowed frequency response 
is now equal to the convolution 

G(u) = F(u) 0 W(u) 

The ideal Fourier transform 
given as: 



of the FIR low pass filter 
F( u ) and W(u) : 

(4.37) 

of a two dimensional LPF is 



F^CujV) = { 



|u^+v^ I <A^ 
I u^+v^ I >A^ 



(4.38) 



The infinite impulse response f 2 (x,y) of this filter has 
to be truncated using a window function W 2 (x,y). The finite 
impulse response filter is given as: 



g2(x,y) = f2(x,y) • 



(4.39) 



The windowed frequency response of this filter is given as: 



G2(u,v) = F2(u,v) 0 



(4.40) 



A good window should have the width of its frequency response 
much smaller than the width of the filter's frequency response 

r2 — 7 

W(u) and should be zero for values of /u +v much 

smaller than A. 
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Near a discontinuity in the filter's frequency response 
F(u), i.e.,' u = G(u) is essentially a convolution of 

W(u) and a step function: 

G(u) = H(u) 0 W(u) (4.41) 

Similarly near a discontinuity in the filter's frequency 
» . 2 2 2 

response F^CujV), i.e., u +v = A , G 2 (u,v) is essentially a 
convolution of W 2 (u,v) and a step function: 

G2(u,v) = H^CujV) 9 W2(u,v) (4.42) 

From Huang's theorm given in (4.33), G 2 (u,v)=G(u) if only 
2 2 

w^Cx^y) = w(/x +y ). 

This means that if w(x) is a good window, then, W 2 (x,y) 
is also a good two dimensional window. 

If other words, a good two dimensional window can be 
derived from a good one dimensional window by extending 
it into a circularly symmetric function. 

C. TWO DIMENSIONAL WINDOW FIGURE'S OF MERIT 

The window's figures of merit found in the one dimen- 
sional case are now applied to the two dimensional case. The 
same figures of merit that may be used for comparison among 
windows, will be re-examined to determine if there should 
be any changes for the two dimensional case. 
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1 . Equivalent Noise Band Width (ENBW) 



The ENBW in one dimensional case is defined as the 
width of a rectangular filter with the same peak power gain 
that passes the same noise power. 

A rectangular filter in two-dimension represents a 
cross-section of a cylindrical filter. The noise power 
passed by the area under the rectangular filter, becomes 
noise power passed by the volume under the cylindrical filter. 

If the peak power gain is the same, then the volume 
divided by the height, results in the area covered by the 
base of the cylinder, which is a circle. 

Two-dimensional ENBW is defined as the circumference 
of the horizontal cross-section of a cylindrical filter, 
that passes the- same noise power. 

The vertical cross-section of a two dimensional 
window is equal to a one dimensional window. Consequently 
the width of the equivalent rectangular window is equal to 

the diameter of the equivalent cylindrical window. Then 

/ 

the two dimensional ENBW is given as: 

(ENBW)2J3 = 7t(ENBW)^j^ (4.42) 

For the two dimensional case, the definition of ENBW 
should be changed from the width of a rectangular to the 
circumference of a cylinder, and the value can be calculated 
directly from the one dimensional case . 
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2 . Coherent Gain 



Two dimensional coherent gain is best defined using 
two dimensional sinusoidal signal. The input signal is 
defined as: 

f(x,y) = A + n(x,y) (4.43) 

^ ? y 

where n(x,y) is two dimensional white noise. 

The output signal or the windowed signal is given as : 

f (x,y) = Aw(x,y)e^ ,y)n(x,y) (4.44) 

w 

The two dimensional frequency response of the win- 
dowed signal is given as: 



F (x,y) 
w 






a a 

= A / / w(x,y)dxdy. (4.45) 

-a -a 

The amplitude of the windowed signal frequency response 
is proportional to the original amplitude A. The propor- 
tionality factor is the double integral of the window function. 

If the window function is given in circularly symme- 
tric form, then the proportionality factor becomes: 



a 


a 


a 


2tt 


/ 


/ w(x,y)dxdy = 


/ 


/ w(r)r d6dr 


a 


-a 


o 


o 



61 



(4.46) 



a 2it 

= / r w(r)dr f d6 

o o 

a 

= 27t f V w(r)dr 
o 

The proportionality factor for a rectangular window 
2 

is equal to ira . The proportionality factor for the other 
window is reduced due to the window smoothly going to zero 
near the boundaries . 

If the proportionality factor in (4.46) is normalized 
2 

by factor ira (that is, divided by the area covered by the 
base of the window) , then the proportionality factor becomes 
coherent gain: 

2 ^ 

COH.GAIN = -y / r w(r)dr (4.47) 

a o 

The coherent gain given in (4.47) above is used as 
window's figure of merit. 

3 . Trade-off Parameters 

Two dimensional trade-off parameters are similar in 
definition with the ones used in the one dimensional case. 

The trade-off parameter a^, a^, b and d are derived from the 
frequency response of a window, which is found using the 
Hankel transform. 

The Hankel transform allows a two dimensional trans- 
form to be found from a one dimensional manipulation. Con- 
sequently the definitions of the trade-off parameters doesn't 
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I 




change. The graphical interpretation of the parameters 
shown- in Fig. 4.1. 



are 



63 



w(r) 





Plot Of Kaiser Window, Used To Illustrate The 
Definition Of Trade-Off Parameters a^ , a2 , b and d 
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V. TWO DIMENSIONAL WINDOW FUNCTIONS 



The two dimensional extension of the previous window 
functions will be given in this chapter. The window's fre- 
quency response is found using Hankel transform. Unfor- 
tunately it is hard to present the Hankel transform of a 
window in a closed form. Only two transforms are presented 
in closed form, they are Rectangular and Barlett windows, 
the rest are presented in their original integral form. 

Without loss of generality, the window function is given 
as w(r) , a two dimensional circularly symmetric function, 
restricted to the interval |r| ^ a, while the window is set 
to unity at the origin. 

00 

w(o) = / f W(f) df = 1 (5.1) 

— 00 

w(r) = o for |r| > a (5.2) 

The frequency response of this window is given as: 

00 

W(f) = / r w(r) J (fr) dr (5.3) 

o ° 

The frequency response consists of a main lobe and side 
lobes which posses a large amplitude ratio in a linear scale. 
The frequency response plot is given in log-frequency scale 
versus log-normalized amplitude (W(f)/W(0)). 
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A. RECTANGULAR WINDOW [2] 



The two dimensional version of the circularly symmetric 
rectangular window is given as follows: 



This window is cylindrical, with the height equal to one and 
the diameter of the base equal to 2a. The vertical cross- 
section of this window is equal to a rectangular window. 

The frequency response of this window is given by the 
Hankel transform as : 



The recurrence formulae is used to solve the integral 
in (5.5a), where: [10] 



1 for I r I _< a 



w(r) = { 



(5.4) 



0 for |r| > a 



a 



W(f)= /rJ(fr)dr 
o 



(5.5a) 




^ {x J (x)} = X J (x) 

IX 1 o 



Then the integral form in (5.5a) reduces to: 



W(f) = ^ 



a 



2 



/ (fr)J (fr) d(fr) 
o 




a 



f 



o 



o 







[ 



4 





(5.5b) 



2 J^(af) 

= a 5 . 

af 

J, (af) 

The window frequency response { — ^ } has a form similar to 
form. The value of this function at the origin, 
f=0, is one. The envelope of the window frequency response 
crosses the frequency axis, whenever J^(af)=0. The first 
zero-crossing occurs at the value of af=3.8. If the value 
of a is increased (meaning the window impulse response is 
lengthened) , the zero-crossing occurs at smaller value of 
f (meaning the window frequency response is narrower) . If 
the value of a is decreased (meaning the window impulse 
response is shortened) , the zero-crossing occurs at larger 
value of f (meaning the window frequency response is wider) . 

At the value of a=0.5, the first zero-crossing occurs at 
the value of f=7.6. 

The ENBW of this window with the value of a is given as 
a=0.5: 



ENBW = .Tr(l) = 3.14 

The COH.GAIN is given as: 

nntj OAT1.T 2 it , 0.78 , p, 

COH.GAIN = 7T f r dr = a no ~ 1*0 

/ r \ ^ U«/0 

TT ( 0 . 5 ) O 

The window and its frequency response are given in 
Figure 5.1(a) and (b) respectively. The main lobe width 
(the width between the origin and the first zero-crossing) is 
7.6. The first side lobe level in 18 dB down from the main 
lobe peak and the side lobe fall-off rate is 9 dB/octabe. 
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-0.5 000 0.5 r 

(a) 




Figure 5.1 

(a) 2-D Rectangular Window, (b) 2-D Spectral Window 
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B. BARTLETT WINDOW 

The two dimensional circularly zymetric Bartlett window 
is given as: 

1 . M 1^1 ^ 

w(r)={ ^ (5.7) 

0 elsewhere 



The shape of this window function is similar to a cone, 
with the height equal to one, and the base in a circle with 
a radius r=a. The vertical cross-section of this window is 
equal to a one-dimensional Bartlett window. 

The frequency response of this window is given by the 
Hankel transform as follows: 

W(f) = / r(l - J (fr) dr 



a a 2 

/ rJ(fr)dr- / — J (fr) dr 

o a. o 

o o 



J (af) ^ a d (fr)J^(fr) 

a - — / r 2 “ 

^ ^ o f 



J-,(af) (af ) , a (fr)J, (fr) 

a - a . 1 / -1 dr 



, a (fr) J- (fr) 

1 / 1 d(fr) 

o f r 



, a d(fr) J^(fr) 

± / ± 

a ^ .4 

o f r 
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i J 2 <af) - ^ f J 2 (fr) d 

o f r 



J (af) , a 

2 — 3 ^ J2(fr) d(fr) 

f af o 



J (af) , a 

-^- 75 — + { / J (fr) d(fr) - 2J, (af)} 

f af o 



af 



^ {af J„(af) - 2 J^(af)} + / J (fr) d(fr) 

O Z J_ JT O 



af o 



af 



1 2 a 

- 7 T {-af J (af)} + — - — ^ / J (fr) d(fr) 

o O / j=- ^ 'J O 



(af ) o 



W(f) = a^{ — ^ J (af) + — ^ / J (fr) d(fr)} (5.8) 

4 0 



(af ) 



(af ) o 



If p=af is substituted into (5.8) then the end result is 



2 r -3 



-2 



W(p ) = a {p / J (t)dt - p J (p)} 
a ^ o o 

o 



(5.9) 



where t is a dummy variable of integration. 

The window and its frequency response is given in Figure 
5.2(a) and (b) respectively, where a=0.5. The ENBW and 
COH.GAIN parameters are given as: 



ENBW = 77(1.33) = 4.18 

0.5 



COH.GAIN = 



2tt 



/ r 



M(0.5) o 






0 . 26 
0 .78 



= 0.34 
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i w(r) 





Figure 5.2 

(a) 2-D Barlett Window, (b) 2-D Spectral Window 
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The main lobe width is about 11.5, The first side lobe 

level is 34 dB down from the main lobe level, and the side 

lobe fall off rate is 12 dB/octave. 

P 

The term / J (t)dt in (5.9) is hardly a closed form, 
c ° 

but it is tabulated. Then the rectangular and relatively 
the Bartlett windows are two-dimensional windows that have 
a closed-form frequency response. 

C. COSINE WINDOW 

The third two-dimensional window to be discussed is a 
cosine window. The two dimensional circularly symmetric 
window is given as: 

cos r for |r| < a 

w(r ) = { ^ (5.10) 

o elsewhere 

The shape of this window is similar to a half sphere and 
the height at r=o is one. The base of the window is a circle 
with radius r=a. 

The frequency response of the window is given as a 
Hankel transform with no closed form. 

^ TT 

w(f) = / r cos r J (fr) dr (5,11) 

2a o 
o 

The effort to derive a closed form is given in the Appendix. 
If a=0.5, then the window and its frequency response are 
given in Figure 5.3(a) and (b) respectively. The parameters 
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rH : w(r ) 
o , 
o 1 




Figure 5.3 

(a) 2-D Version Of Cosine Window (cross-section) , 
(b) 2-D Spectral Window 
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ENBW and COH.GAIN are given as: 



ENBW = tt(1.33) = 3.86 

COH.GAIN = — — TT f r cos ttp dr = = 0-'+6 

/ r\ n \ ^ U*/0 

tt( 0 . 5) o 

The main lobe width is 10.5 while the first side lobe 
level is 26 dB down from the main lobe and the side lobe 
fall off rate is 12 dB/octave. 

D. HANN WINDOW 

The two dimensional circularly symmetric Hann window 
is given as: 

cos^ -^r = 0.5+0. 5 cos — r |r| <a 

w(r) = { ^ (5.12) 

o elsewhere 

The vertical cross-section of the window is equal to a 
one dimensional Hann window. 

The frequency response is given by the Hankel transform 

as : 

J (fr) dr 
o 

cos — r J (fr) dr (5.13) 
a o 



o 

w(f) = / r(0.5+0.5 cos — r) 

3 . 

a 

J (af) a 

= 0.5 a — p +0.5 f 

o 
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■m 



£ 




The first term in (5.13) has a closed form as in the case 
of a Rectangular window. The second term has no closed form 
as in the case of the cosine window. 

For a window length of a=0.5, the window and its frequency 
response are given in Figure 5.4(a) and (b) respectively. 

The ENBW and COH.GAIN parameters are given as: 

ENBW = IT (1.5) = 4.7 

COH.GAIN = f r cos^ Trr dr = = 0.2 9 

/ r\ c \ ^ U#/o 

-n-( 0 . 5) o 

The main lobe width is 13.7 and the first side lobe level 
is 34 dB down from the main lobe level. The fall-off rate 
of the side lobes is 18 dB/octave. 

E . HAMMING WINDOW 

The two dimensional circularly symmetric Hamming window 
is given as; 

0.54 + 0.46 cos — r |r| < a 

w(r) = { ^ (5.14) 

o elsewhere 

The shape of this window is a smooth cone sitting on a 
cylindrical base. The vertical cross-section of the window 
is equal to the one dimensional case. 

The frequency spectrum of the window is given by the 
Hankel transform as: 
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I 





20 log^g|W(f)/W(0) I 

b=13.0 50. log f 







-100. _ (b) 

Figure 5.4 

(a) Von Hann Window in 2-D, (b) Spectral Window 
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w(f ) 



/ r(0.54 + 0.46 cos — r) J (fr) dr 
o 



J., (af ) a 

0.54 a 5 + 0.46 / r cos — r J (fr) dr 

r ^ a o 

o 



(5.15) 



The first term in (5.15) is the familiar frequency re- 
sponse of the rectangular window and the second term is the 
response of a cosine window. 

If the length of the window is a=0.5, then the window 
and its frequency response are given in Figure 5.5(a) and (b) 
respectively. The ENBW and COH.GAIN parameters are given as: 



ENBW = tt(1.36) = 4.3 

COH.GAIN = — — ^ / r(0.54 + 0.46 cos ■;n^)dr=§^4l=0 • 36 

IT ( 0 . 5 ) o 

The main lobe width is equal to 13.5 and the first side 
lobe level is 47 dB down from the main lobe level. The side 
lobe fall off rate is 6 dB/octave. 

F. BLACKMAN WINDOW 

The two dimensional circularly symmetric Blackman window 
is given as: 

0.42 + 0.5 cos-^ + 0.08 cos— r |r| < a 

w(r) = { ^ ^ (5.16) 

o elsewhere 
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Figure 5.5 

(a) Hamming Window in 2-D, (b) Spectral Window 
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The vertical cross-section of this window is equal to 
the one-dimensional Blackman window. 

The window frequency response is given by the Hankel 
transform: 



W(f) = 



/ r {0.42 + 0. 5 cos- r + 0.08 cos — r}J (fr)dr 
o ^ a o 



J (af) a 

= 0.4 2 a 5 +0.5 / r cos-^ J (fr)dr 

f a o 

o 



+ 0.08 



a 

/ r cos 
o 



— r J (fr)dr 
a o 



(5.17) 



If the extent of the window given as a=0.5, the window 
and its frequency response are given in Figure 5.6(a) and (b) 
respectively. The ENBW and COH.GAIN parameters are given as: 



ENBW = 7t(1.57) = 4.9 



COH.GAIN = 



2ir 



(0.5) ir o 



0.5 

/ r (0.42 + 0.5 coSn - V -+0 . 08 cos 7 ^-^r)dr 
I U . o U • o 



17 

0.78 



= 0.22 



The main lobe width is equal 18.5, which is almost triple 
of that of the rectangular window. The first side lobe level 
is 60 dB down from the main lobe and the side lobe fall off 
rate is 9 dB/octave. 
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Figure 5 . 6 

(a) 2-D Blackman Window, (b) Spectral Window 
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6. KAISER WINDOW 



Two dimensional circularly symmetric Kaiser window is 
given as: 



I [S 

o a 

rm 

o 



w(r) = { 



r < a 



elsewhere 



(5.18) 



The vertical cross-section of the window is equal to a 
one dimensional Kaiser window. 

The frequency spectrum of the window is given by the 
Hankel transform: 



a I [6 

W(f) = / r -2 — 

o 






J (fr)dr 
o 



(5.19) 



which cannot be reduced to a closed form. 

If a=0.5 and 6=6 then the window and its frequency res- 
ponse are given in Figure 5.7(a) and (b). The ENBW and 
COH.GAIN parameter are given as: 

ENBW = ir(1.5) = 4.71 



COH.GAIN 



2tt 

tt( 0 . 5) 



0.5 

/ 

o 



rm 



dr 



0.48 

0.78 



0.61 



The main lobe width is 14.0 and the first side lobe level 
is 44 dB down from the main lobe peak. The side lobe fall off 
rate is 6 dB/octave. 

In this case, the window's characteristics depend strongly 
on the values of a and 6. With 6=6, the window is similar to 
a Hamming window . 
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(a) 2-D Kaiser Window, (b) Spectral Window 
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TABLE VI. 1 



Two-Dimensional Window Functions and Figures of Merit 





WINDOW 


ENBW 


COH.GAIN 


^1 


^2 


b 


d 










dB 


dB 


Hz 


dB/octave 


1. 


Retangular 


3.14 


1.0 


-18 


-38 


6.6 


-9 


2. 


Bartlett 


4.18 


.34 


-34 


-58 


11.0 


-12 


3. 


Cosine 


3.86 


.46 


-26 


-55 


10.0 


-12 


4 . 


Hann 


4.71 


.29 


-34 


-64 


13.0 


-18 


5. 


Hamming 


4.3 


.36 


-47 


-51 


13.0 


-6 


6. 


Blackman 


4.9 


.22 


-60 


-66 


19.0 


-9 


7. 


Kaiser 


4.7 


.61 


-44 


-61 


14.0 


-6 
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VI. CONCLUSION 



One of the most straightforward applications of the rec- 
tangular window is to truncate an infinite length impulse 
response. This can be thought as multiplying the impulse 
response by a window function that is unity in some finite 
interval between ±a, and zero elsewhere. This multiplica- 
tion is equivalent to convolving the Fourier transform of 

s inx 

the impulse response with a function of the form — ~ — . This 
causes ringing in the original frequency response. Increasing 
the length of the finite impulse response (by increasing the 
length of the rectangular window) decreases the duration of 
the ringing, but has no effect on the peak overshoot. This 
effect is called the Gibbs phenomenon. 

A considerable amount of effort has been put into finding 
a good window function, i.e., window functions which reduce 
the peak overshoot of the Gibbs oscillation. It is generally 
agreed that the frequency response of a good window should 
have a narrow central lobe (with a large amplitude) and side 
lobes with comparatively small amplitude. There are many 
window functions readily available in the literature, however 
most window functions are given only in one dimension. The 
one dimensional window functions and their Fourier transform 
are given in detail in Chapter III of this thesis. Figures 
of merit for one dimensional window are powerful measure and 
are used to compare the characteristic of the various windows. 
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It has been the objective of this thesis to examine two 
dimensional window functions. The two dimensional window 
functions presented here have been a direct extension of 
corresponding one dimensional window functions. The two 
dimensional window function is formed from a one dimensional 
window function, using the relation W 2 (x,y)=w(/x +y ). This 
obviously rotates the one dimensional function about the 
origin in the x-y plane. However it does not simply rotate 
the frequency response of the one dimensional about the origin 
in the u,v plane. The two dimensional frequency response of 
the window can be calculated using the Hankel transform. The 
principal advantage of this approach is that the frequency 
transform of two dimensional window can be found using a one 
dimensional manipulation. 

The two dimensional window figures of merit are a direct 
extension of the one dimensional case. However the definition 
of the ENBW is changed from the width of a rectangular filter 
in the one dimensional case to the circumference of a circular 
filter in the two dimensional case. Table VI . 1 shows the 
comparison between one dimensional and two dimensional window 
figure of merit. From the table it is clear that the two 
ENBW parameters are related by the equation (4.42). 

The coherent gain decreases slightly except for the Kaiser 
window where it is higher, and the rectangular window where it 
is constant. The first side lobe level, which is shown by the 
parameter a^ , is lower for the two dimensional windows, however 
the width of the main lobe is much greater in each case. The 
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side lobe fall-off rate stays the same except for the Blackman 
window where it is slower and the rectangular window where it 
is faster. 

In conclusion, while one dimensional and two dimensional 
windows can be simply related to one another in the time 
domain, no direct, general parallels exist to relate their 
transforms in the frequency domain. 
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TABLE VI. 1 



Comparison Between Two Dimensional And 
One Dimensional Window Figures of Merit 





ENBW 


COH. 


GAIN 


a 


1 




b 


d 




WXi>i JJU w 


2D 


ID 


2D 


ID 


2D 


ID 


2D 


ID 


2D 


ID 


Rectangular 


3.14 


1.0 


1.0 


1.0 


-18 


-13 


6.6 


0.81 


-9 


-6 


Bartlett 


4.18 


1.33 


.34 


.5 


.34 


-26 


11.0 


1.63 


-12 


-12 


Cosine 


3 . 86 


1.23 


.46 


0 . 64 


-26 


-23 


10.0 


1.35 


-12 


-12 


Hann 


4.71 


1.5 


.29 


0.5 


-34 


-32 


13.0 


1.87 


-18 


-18 


Hamming 


4.3 


1.36 


.36 


0 . 54 


-47 


-43 


13.0 


1.91 


-6 


-6 


Blackman 


4.9 


1.57 


. 22 


.46 


-60 


-58 


19 . 0 


2.82 


-9 


-18 


Kaiser 


4.7 


1.5 


.61 


.49 


-44 


-43 


14.0 


2.2 


-6 


-6 
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APPENDIX 



The effort to put a Hankel transform of a cosine function 
in a closed form is given below. The cosine function arises 
from Hankel transform of a cosine window. 

The Hankel transform of a cosine function is given below: 

^ 7T 

w(f) = / r cos r J (rf) dr Irl < a (A.l) 

2a o ' ' — 

o 



The above integral form has no known closed form, as proved 
by the following discussion. 

Three different methods are used to bring the integral 
into a closed form, they are: straight integration, partial 

differential and table of integration. 



A. STRAIGHT INTEGRATION 

To use straight integration method to solve the integra- 
tion, the variable needs to be changed into trigonometric 
variables . 



Let : 



9 = -1^ 



TTr 

2a 




6 . 



• 7T 

The interval of the integral becomes, 9=o and for r=o 

and r=a respectively. 

Substitution of the new variables into Hankel transform 
gives : 

, 2 it/ 2 

w(f) = (^) / e cos 9 J (X9) d9 (A. 2) 

f o 

o 
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where X 



= 2(^). This new integral form gives more pleasing 
appearance than the previous one. 

Substitute the Bessel function J^(X0) inside the integral 
with Bessel equation, given as follows: 



J "(x) + - J ' (x) + J (x) = 0 
O X o o 



(A. 3) 



Let X = X6, and do some manipulation. The Bessel equation 
becomes : 



J (X0) = - i {J "(A0) + J- J 'U0)} 
o o 6 o 



(A. 4) 



The integral form in A. 2 becomes: 



u/ 2 



w(X) = — / 0 cos 0 {J "(X0) + ^ J '(X0)}d0 

f2 o o 0 o 



1 ir/2 -n/2 

— s-{0 cos 0 J '(X0) I / {(cos 0-0 sin 0)J '(X0) 
-P ^ o o 

t o o 



- cos 0 J ' (X0)}d0} 
o 



ir/ 2 






f 0 


sin 0 J * ( X0 ) 

/-V 


d0} 


o 


CJ 






tt/ 2 


tt/2 


sin 


CD 

) 

>" 

CD 

1 


/ 




o 


o 



1 rU , /TT 



tt/2 



(tt J (^-X)} + —TT f sin0J(X0)d0+ — ^ w(X) 



'71^2 71 



f o 
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(1 - ^)w(X) 



2t 




ir/2 

7f / sin 9 J (X9) d9 
^ o 



The end result is: 



w(X) = -V^ / sin 6 J (X9) d9 



x2_i ir/2 




J (5-X) (A. 5) 

O 2 



By replacing J^(X9) inside the integral as before, it may 
be continued at will, but the term with the integral form 
cannot be eliminated. The conclusion is, the integral form 
in A. 2 cannot be put in a closed form by this method. 

B. PARTIAL DIFFERENTIAL METHOD 

Return to the integral form given in A. 2: 



X 9 ' 

W(X) = (^r) f 0 cos 9 J (X9) d9 



Define new equations: C(X) and S(X) 



it/ 2 

CtX) = / 0 cos 9 J (X9) d9 

o 



(A. 6) 



it/ 2 

S(X) = / 0 sin 9 J (X9) d0 

o 



(A. 7) 



Differentiate to X as C'(X) and S’(X) 



C’ (X) = ^ C(X) 



tt/2 „ 

f 9 cos 9 J ' (X0) d9 
o 



(A. 8) 
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(A, 9) 



S’ (A) 



= ^S(A) 



tt/2 ^ 

/ 0 sin 0 J ’(A0) d0 

o 
o 



Multiply both A. 8 and A. 9 with A: 



ir/2 ^ 

AC’(A) = / 0^ cos 0 J ’(A0) d(A0) 

o 

o 



(A. 10) 



tt/2 r . 

AS’(A) = / 0 sin 0 J '(A0) d(A0) (A. 11) 

o 

o 



Using Bessel identity ’ ( A0 )d( A0 ) =dJ^( A0 ) , substitute 
it into A. 10 and A. 11: 



tt/2 „ 

AC’(A) = / 0^ cos 0 dJ (A0) (A. 12) 

o 

o 



AS' (A) 



tt/2 r . 

f 0 cos 0 dJ (A0) 
o 



(A. 13) 



Expand both integral in A. 12 and A. 13 into partial integration: 

„ ir/2 tt/2 ^ 

AC'U) = 0 cos 0 J (A0) I - / (20COS0-0 sin0)J (A0)d0 

o ' _ o 



AS' (A) 



tt/2 ^ 

- / (20 COS0-0 sin0) J (A0) d0 (A. 14) 

o 

o 



0 sxn0 J (A0) (20 sin0+0 cos0)J (A0)d0 

o ' o 

o 
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(A. 15) 



= ^ J (Ja) - / (20 sin0+9^ cos0)J (Ae)d0 

2 o 2 o 

o 

Rewrite both equation in A. 14 and A. 15, using definition 
of C and S. It will end up with the following pair of differ- 
ential equations: 

AC’ (A) = -2C(A) + S’ (A) (A. 16) 



AS’ (A) = ?• J (?-A) - 2S - C’ 
I O I 



(A. 17) 



Rearrange the last two equations: 

AC’ + 2C + S’ = 0 (A. 18) 

AS’ + 2S + C’ = 5- J (Ja) (A. 19) 

The last two pairs of equations show that it may give a single 
equation in terms of S or C. And that means it will end up in 
integral form of C or S. The conclusion is, the integral 
form in A. 2 cannot be put in a closed form by this method. 

C. TABLE OF INTEGRATION 

Taylor series expands any analytic function f(x) into a 
power series of x. Consequently cos x may be expanded into 
a power series of x : 
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2 

cos X = f(x ) 



(A. 20) 



Substitute A. 20 into A. 2; 



w(x) = / X f(x ) J (fx)dx 

o ° 



(A. 21) 



where we change the integration variables with the dummy 
variables x and f. 

From Table of Integration [Ref. 10], the integral form 
shown below: 



(x) -ft J (t)dt 

2k o 

o 



(A. 22) 



can be written in a closed form. The term t(i (t) in the 

o 

above equation is replaced using Bessel equation: 

tJ (t) = -tJ "(t)-J ’(t). 
o o o 

The integral form in A. 22 becomes: 



(x) = / t^^{-tJ "(t)-J '(t)} dt 

2k ^ o o 

o 



/ '(t)dt- / t^^J '(t)dt 

o o 

o o 



-x^^''’^J »(x) + (2k+l) / t^^J '(t)dt- f t^^J '(t)dt 
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_ _^2k+lj i(x)+2k{x^^J (x)-2k / (t)dt} 

o o o 

o 

I.,(x) = 2kx^^J (x)-x^^'^^J '(x)-4k^I„, ,(x) (A. 23) 

ZK O O Zs^"^ jL 



The last equation shows that l 2 j^(x) reduces to lower order 
^2k This form is called reduction form. For k=o, 

A. 22 and A. 23 become: 



I (x) = / t J (t)dt = =xJ ' (x) 
o o o 

o 



(A. 24) 



The equation in A. 24 agrees with the property of Bessel 
function . 

If the reduction form applied repeatedly and written in 
a general form, A. 23 may be written as: 



2 k. 



X 

/ t"^tJ (t)dt = 2k.Poly(2k) . J (x)+x.Poly(2k) . J ’(x) (A. 25) 

o o o 

o 



Lets return to equation A. 20 



2 2k 

cos X = f(x ) = Z ^k ^ 

k=o 



(A. 26) 



Substitute A. 26 into A. 25: 



/ f(t )t J (t)dt 
o 



Z a. x^^ . J (x)+ Z b, ' (x) 

k=o ^ ° k=o ^ 



(A. 27) 
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Apply differentiation to both left and right hand side of A. 27 



I f^x^ (x) = J (x) Z 2ka,x^ ^+J ’ (x) Z (2k+l)b 

1 k o o 1 _ k o , _ 

k=o k=o k=o 



+ J ’ (x) Z a. x^^ + J "(x) Z b, x^^^^ 
o k o k 



Replace J^"(x) 

J "(x) = - 
o 



using Bessel equation: 

x"^J ' (k) - J (x) 
o o 



(A. 28) 



Z f (x) = J (x) { Z 2ka x^^"^ - Z b, 

k,o X o o ^ k k 



+ J '(x) {Za,x^^ + Z (2k+l)b, x^^-Zb, x^^} 

O X XX 



J (x) { Z (2k+l)a, X - Z b,x 
o kt k 

o o 



+ J ’ (x) {Z (a +2kb, )x^^} 
o o ^ ^ 



Equate the coeficient of the equal term and ’ from the 
left hand side and right hand side: 



0 = E (a, +2kb, )x^^ 
k k 



Therefore : 



^k-l^^k-1 

f, = (2k+l)a, , -b ->- a = 7 ^ 

k k+1 k k 2k 



0 = 



' ■ 2k 



These last two equations give the solution to the 
of the cos X form. Lets expand cos x into power 



cos X = f 



+ f ^x 



f 2^ 



" ^ 3 " 



then using relations given in (A. 30) and (A. 31) 



a =0 
o 



b = 0 



b = - 






f,-f /4 
1 o 



3 

etc 



32 



f,-f /4 
1 o 

16 

f,-f /4 
1 o 

192 



(A. 29) 
(A. 30) 



(A. 31) 

(A. 32) 

integral 
series : 
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The conclusion is, we succeeded in solving Hankel trans- 
form of cos X function, but we failed to express it in a 
closed form. 
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